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We introduce a new class of unitary transformations based on the su(l, 1) Lie algebra that genera- 
lizes, for certain particular representations of its generators, well-known squeezing transformations 
in quantum optics. To illustrate our results, we focus on the two-mode bosonic representation and 
show how the parametric amplifier model can be modified in order to generate such a generalized 
squeezing operator. Furthermore, we obtain a general expression for the bipartite Wigner function 
which allows us to identify two distinct sources of entanglement, here labelled by dynamical and 
kinematical entanglement. We also establish a quantitative estimate of entanglement for bipartite 
systems through some basic definitions of entropy functionals in continuous phase-space represen- 
tations. 

PACS numbers: 02.20.Qs, 03.65. Ca, 03.67.Mn 



I. INTRODUCTION 

In the last few years, physics has experienced the ap- 
pearance of two relatively young branches with strong 
appeal in both theoretical and experimental aspects. La- 
belled by Quantum Information and Quantum Computa- 
tion, these branches apparently entangled have attracted 
since then a lot of attention from researchers working 
in well-established areas in physics (such as, for ins- 
tance, solid state physics, nuclear physics, high energy 
physics, general relativity, and cosmology [l|). The in- 
terdisciplinarity provided by Quantum Information and 
Quantum Computation is basically focussed upon fun- 
damental physical concepts that constitute the corner- 
stones of quantum mechanics. Hence, concepts related 
to nonclassical states, superposition principle, entangle- 
ment, phase-space representations, quantum teleporta- 
tion, quantum key distribution, among others, represent 
nowadays common words in many scientific papers cove- 
ring different areas of knowledge in physics. 

In particular, let us restrict our attention to quan- 
tum information theory and its description by quantum 
continuous variables [2j, where entanglement effects can 
be efficiently produced in laboratory through the ade- 
quate manipulation of continuous quadrature amplitudes 
of the quantized electromagnetic field. In this promi- 
sing scenario, the squeezed light [H, 0] has a prominent 
role in the experimental implementation of continuous- 
variable entanglement, since the degree of imperfection in 
entanglement-based quantum protocols depends on 'the 
amount of squeezing of the laser light involved'. From a 
conceptual point of view, some questions related to the 
entanglement measures (or separability criteria) still re- 
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main open in the specialized literature Q. For instance, 
entanglement of formation 'ja\ and concurrence Q are 
now widely accepted as entanglement measures for the 
two-qubit case. However, it is worth noticing that diffe- 
rent approaches to this problem exist which are outlined 
by means of information-theoretic arguments Q . 

The main goal of this paper is to present some contri- 
butions to certain specific topics in quantum information 
theory that allow us to go further in our comprehen- 
sion on the entanglement process in ideal bipartite sys- 
tems (the interaction with any dissipative environment 
@ is discarded in a first moment). For this purpose, 
we first construct a general family of unitary transfor- 
mations associated with the su(l, 1) Lie algebra genera- 
tors that generalizes - if one considers the one- and two- 
mode bosonic representations of its generators - two well- 
known expressions of squeezing operators Q . Following, 
we study a modified version of the parametric ampli- 
fier (lol | with emphasis in obtaining the solutions of the 
Heisenberg equations and its respective time-evolution 
operator. In particular, we show the efficacy of this model 
in generating the generalized two-mode squeezing opera- 
tor through its connection with the time-evolution ope- 
rator. 

The next step then consists in performing a prelimi- 
nary study via Wigner function on the qualitative as- 
pects of entanglement for a general class of bipartite sys- 
tems whose dynamics is governed by arbitrary quadratic 
Hamiltonians. As a by-product of this investigation, we 
obtain a general integral representation for the bipartite 
Wigner function that leads us to identify two distinct 
sources of entanglement (here labelled by dynamical and 
kinematical entanglement). The last contribution refers 
to a direct application of the results obtained by Piatek 
and Leohski [Lj| for the intermode correlations in con- 
tinuous phase-space representations, which are based on 
the Wehrl's approach [12| regarding some definitions of 
entropy functionals and their inherent properties. In fact, 
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we introduce a correlation functional that permits us to 
measure the degree of entanglement between both parts 
of the joint system for any initial conditions associated 
with the Hamiltonian operator. It is important to em- 
phasize that the sequence of topics covered in this work 
presents an inherent logical consistency that improves our 
comprehension on the subtle mechanisms associated with 
the entanglement effects in bipartite systems. 

This paper is structured as follows. In Section II, we 
present a mathematical statement that leads us to cons- 
truct a general family of unitary transformations associ- 
ated with the su(l, 1) Lie algebra where, in particular, 
the two-mode bosonic representation is emphasized. In 
Section III, we obtain the solutions of the Heisenberg 
equations for the parametric amplifier model, and show 
the subtle link between time-evolution operator and ge- 
neralized two-mode squeezing operator. The results ob- 
tained are then applied in Section IV, within the con- 
text of Wigner functions, in order to establish an initial 
discussion on entanglement for certain groups of bipar- 
tite states of the electromagnetic field. Section V is de- 
voted to establish a reasonable measure of entanglement 
which is based on some specific information-theoretic ar- 
guments. Besides, we illustrate our results through two 
different examples of initially uncoupled bipartite states 
for the model under investigation. Finally, Section VI 
contains our summary and conclusions. 



II. SQUEEZING OPERATORS ASSOCIATED 
WITH THE sufl, 1) LIE ALGEBRA 



Let us initially introduce the generators of the su(l, 1) 
Lie algebra K± and Ko, which satisfy the following com- 
mutation relations: 

[K_,K+] = 2K and [K ,K±] = ±K± . 

The Casimir operator is defined within this context 
through the mathematical identity 



K 2 



K 2 



(1/2){K_,K + } 



with [K 2 ,K±] = [K 2 ,K ] = 0, where {K_,K+} repre- 
sents the anticommutation relation between the opera- 
tors K_ and K + . Furthermore, let us also consider the 
abstract operator 



deal with some applications of the composition formula 
with emphasis on unitary transformations which resem- 
ble, for certain representations of the generators K± and 
Ko, the squeezing transformations in quantum optics. 

Lemma II. 1 Let T(f2-t,f)o) and T(A±, Ao) be two abs- 
tract operators whose functional forms obey equation 
(QJ). For a given set of arbitrary c-number parameters 
{fi±, Oo, A±, Ao} it is always possible to verify the gene- 
ral composition formula 



T(fi±,n )T(A ± ,Ao) = T(£±,£ ), 



(2) 



where £± and So are solutions of the coupled set of non- 
linear equations 

A+ + B B + [Aq-A+ (A_+5_)] _ £+ 

B (A_+B_) "£_' [6a) 



VBo/Ao (A- + BJ) = (S-//3) sinh(/3), 



(3b) 



with [3 = [(£o/2) 2 — E+E-] 1 / 2 . The c-number functions 
A±, Ao, B±, and Bo, present in the Ihs of equations \3a\) 
and \3b\) . are connected with {£l±, Qo, A±, Ao} through 
the identities U2 l: 



A+ = 



(0±/<£) sinh(» 



coshO) - (n /2<£) sinh(^) ' 
Ao = [cosh(0) - (n /2<£)sinh(</>)]~ 2 . 
(A±/0) sinh((9) 



B 



cosh(<9) + (Ao/20) sinh(0)' 



(4a) 
(4b) 
(4c) 
(4d) 



for 



Bo = [cosh(0) + (Aq/20) sinh(<9)] 2 , 

= [(O /2) 2 -O+O-] 1/2 , 
= [(A /2) 2 - A + A_] x / 2 . 



Proof II. 1 Firstly, we apply the generalized normal- and 
antinormal- order decomposition formulas established by 
Ban flA l for exponential functions of the generators of 
su(l,l) Lie algebra to the abstract operators T(f2-|-,fio) 
and T(A±, Ao), that is, 



T(Q±,n ) = e ^+K +e MA )Ko e A_K_ 



(5) 



and 



T(n±,fi ) :=e i2 +"+ 



_ „f2 + K + +n K +n„K_ 



(1) 



written in terms of the arbitrary c-number parameters 
Cl± and Oqj whose generalized normal- and antinormal- 
order decomposition formulas have already been establi- 
shed in literature [H, [3| • These preliminary considera- 
tions on the su(l, 1) Lie algebra lead us to demonstrate 
an important composition formula involving the product 
of two abstract operators, each one being defined as in 
equation ([T]) and characterized by a particular set of ar- 
bitrary c-number parameters. In fact, this section will 



T(A±,A C 



e B_K_ e ln(B )Ko e B+K + 



(6) 



where the c-number functions A±, B±, Ao, and Bo were 
defined by means of the identities J4a| ]- f^rf[ ). In this way, 
the product T(f2±, fio)T(A±, Aq) can be expressed as 



T(ft±,r!o)T(A ± ,Ao) 



eJ 4 + K +e ln(A )K 0e ( J 4_+B_)K_ 
x e ln(B )K 0e B + K + 
e A + K +e ln(A )Ko e C + K + 
x e hx(Co)K 0e C_K 



(7) 
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with C± and Co given by FZJ /: 



C+ = 

C- = 

c a = 



1 -B B+ (A_ +B_Y 
B (A_+B_) 



1 - B B+ (A. 

B 



[1-B B + (A_+B_)] 2 ' 
The second step consists in using the relationship 

e ln(A )Ko e C+K + = e AoC + K +e ln(^o)K 

for the second and third exponentials in the second equa- 
lity on the rhs of equation with the aim of establishing 
the intermediate result 

T(ft±,r!o)T(A ± ,Ao) = e (A + +AoC + )K +e MA Co)K 0e C_K__ 

The rhs of this equation represents the normal-order de- 
composition of an abstract operator defined as 

T(E±,E ) := e E + K + +E ° Ko+E - K -, 

where the c-number parameters £± and S satisfy the 
following mathematical relations: 



A++A C+ = 



(S+//3) sinh(/3) 



cosh(^) - (So/2/3) sinh(/3) ' 

(S_//3)sinh(/3) 
cosh(/3) - (So/2/3) smh(/3)' 



A Q C'o = [cosh(/3)-(S /2/3)sinh(/3)] 2 , 

with (3 — [(So/2) 2 — S + S_] 1//2 . Consequently, substi- 
tuting the definitions of C± and Co in these relations, 
we obtain a coupled set of nonlinear equations that per- 
mits us not only to establish a link between {S±, So} and 
{Q±, Oo, A±, Ao}, but also to verify the general composi- 
tion formula @). H 

An interesting consequence from this mathematical 
statement is associated with the construction process of a 
general family of unitary transformations where the abs- 
tract operator ([T]) has a central role. To carry out this 
task let us first establish a corollary which is directly re- 
lated to the Lemma II. 1. 

Corollary II.l For Q, + — £, f2_ = — £*, and f2 = iw, 

with £ G C and lu G K, £/ie abstract operator (QJ) repre- 
sents a generator of unitary transformations associated 
with the su(l, 1) Lie algebra. 

Proof II. 2 Basically, the idea is obtaining a specific sub- 
set of arbitrary c-number parameters such that 

T(n±, n ) [T(n±, fi )] f = [T(n±, n rf t(o±, n ) = 1. 

Hence, let us initially investigate under what circumstan- 
ces the mathematical relation 



T(n±,n )[T(n ±> no)] 



is verified. In fact, this condition can be promptly derived 
from the general composition formula 0) for A± = fit-, 
Ao = Qq, and S± = So = 0, with the additional res- 
trictions A± = — f2± and Ao = — f^o ■ Consequently, the 
equalities Q± = — Qt. and Qq = ~^o are satisfied only 
for 17+ = £, Q- = — £*, and f2o = with ( £ C and 
TTiis same particular subset of arbitrary c-number 
parameters can also be obtained from the analysis of 

[T(r> ± ,r» )] t T(r! ± ,r!o) = i, 

which implies that T(£, a;) gives a family of unitary trans- 
formations for a general class of representations associ- 
ated with the su(l, 1) Lie algebra. ■ 

To illustrate our results, let us consider as a first 
example the single-mode bosonic representation of the 
Heisenberg-Weyl algebra where the generators K± and 
K are expressed as K + = (l/2)a^ 2 , K_ = (l/2)a 2 , and 
Ko = (l/2)(ata+l/2), with a and being, respectively, 
the boson annihilation and creation operators satisfying 
the well-known commutation relation [a, a^] = 1. In this 
case, the unitary operator T(£, ui) assumes the form 

Ti(£,w) = e |[^ 2 +i W (ata+l/2)-£*a 2 ] ) (g) 

which coincides with the squeezing operator 
S(0 := e K? at2 -C*a 2 ) whcn w = o 

It is worth mentioning that, in particular, equation ([5]) 
can be promptly used to derive the quantum analogous 
of the Fresnel transform in classical optics [l5l |. 

Another interesting example encompasses the two- 
mode bosonic representation in which the generators 
are specifically given by K + = a^b^, K_ = ab, and 
K = (l/2)(at a + btb+ 1), where [a, at] = [b,b+] = 1. 
In this context, the unitary operator 

T 2 (C,U>) = e Ca t bt+i( w /2)(at a +btb+l)-C*ab /qn 

recovers the two-mode squeezing operator 

S 2 (0 = exp(^a t b t - £*ab) for w = 0. 

Moreover, the action of T 2 (£, ui) in the annihilation ope- 
rators for each mode of the electromagnetic field implies 
in the following results: 

T 2 (£,w)aT 2 (£,cj) = [cosh(0) + i(w/2<£) sinh(0)] a 

+ (C/0)sinh( ( /»)bt, (10) 

T 2 (e,w)bT 2 (C,w) = [ooflh(^) + i(w/20) sinh^)] b 

+ (£/0)sinh(</>)a t , (11) 

being = [|£| 2 — (w/2) 2 ] 1 / 2 . This set of unitary transfor- 
mations for bj — has its counterpart in the quantum des- 
cription of physical processes involving parametric ampli- 
fication [HI [li|. Recently, Pielawa et al. [17] have pro- 
posed a new method for generating two-mode squeezing 
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in high-Q resonators using a beam of atoms (which acts 
as a reservoir for the field) with random arrival times. In 
particular, the authors have used the unitary two-mode 
squeezing operator S2(£) to bring an effective Hamilto- 
nian - which describes resonant single-photon processes 
- to the well-known Jaynes-Cummings form, where the 
new bosonic operators are connected to the old ones by 
two- mode squeezing transformations. Another interes- 
ting application is based on the analogy between phonons 
in an axially time-dependent ion trap and quantum fields 
in an expanding/contracting universe, where the multi- 
mode squeezing operator represents the basic mechanism 
for cosmological particle creation . 

It is worth noticing that there are some textbooks on 
quantum optics Q which discuss particular cases (if one 
compares with those exposed here) of one- and two-mode 
squeezing operators in different contexts in physics and 
their connections with nonclassical states of the electro- 
magnetic field. However, the dynamical origin of the real 
parameter uj present in equation @, for example, has 
not been investigated up until now in literature, and this 
fact will be our object of study in the next section. For 
this intent, we obtain the exact solutions of the Heisen- 
berg equations for a specific nonresonant system like the 
parametric amplifier, and show the link between time- 
evolution operator and T2(£,u>). 



III. PARAMETRIC AMPLIFIER 

The parametric amplifier model proposed by Louiscll 
et al. [10J consists basically of two coupled modes of the 
electromagnetic field, which play a symmetrical role in 
the amplification process [IB]. Such dynamical elements 
are usually described by the Hamiltonian (h = 1) 

H(t) =w a a t a + tj b b t b + Kabe i ' )t + K*a t b t e~ i ' ?t , (12) 

where a (b) is the annihilation operator for the signal 
(idler) mode, ry provides the frequency of the pump field 
(which has been assumed strong enough to be expressed 
in classical terms), and k represents a complex coupling 
constant being proportional to the second-order suscep- 
tibility of the nonlinear medium and to the amplitude of 
the pump. Moreover, let us introduce a small deviation 
6 in the usual definition of r/ such that rj = u) a + LUb + 5 
with 8/{uj a + u>b) <C 1 (by hypothesis we are assuming 
that 8 comes from a non-perfect match between the fre- 
quencies 1] and uj a +u>b). In this case, the solutions of the 
Heisenberg equations are given by 

a(t) = e- i ^ +5 / 2 ) t {[cosh(^) + i(V2y) sinh(^)] a(0) 

-i(rc*/V) sinh(^)bt(0)}, 
b(t) = e- i{uJb+s/2)t {[cosh(<pt) +i(6/2tp)smh(<pt)]b(0) 

-i(K*/tp) sinh(<pt)a t (0)}, 



plus their Hermitian conjugates for tp — [\k\ 2 — (S/2) 2 ] 1 / 2 
fixed. These solutions lead us to verify that [a(i), at (£)] = 
[b(i),bt(£)] = 1, which implies in the unitariety of the 
time-evolution operator U(i). This fact leads us to show 
that n a (t) — n&(i) = n a (0) — n&(0) (conservation law) for 
n c (t) := c'(t)c(t) is easily seen to hold; therefore, the 
intensity correlation function (n a (i)n(,(t)) can be written 
as (n a (t)n b {t)) = (n 2 a {t)) + (n a (t) [n 6 (0) - n a (0)]>. For 
instance, when the initial state coincides with the num- 
ber states {\n ai nb)}n a ,n b eNi the second term on the rhs of 
this equation is reduced to (n& — n a )(n a , rib\n a (t)\n a , rib) 
for n a ^ rib. So, if one considers n a — rib, there is no con- 
tribution from this term and consequently, (n a (t)iib(t)) 
corresponds to the maximum violation of the Cauchy- 
Schwarz inequality (a^abtb) < (at 2 a 2 ) for the parame- 
tric amplifier Q. 

Next, we show how the time-evolution operator U(i) 
can be connected with the unitary operator T2(£, u>). To 
carry out this task, let us initially mention that the time- 
dependent global phase factors present in the solutions of 
the Heisenberg equations are obtained through the action 
of the rotation operator 

RK, LUb, 5; t) = e -i*[K+<5/2)ata+K+*/2)btb] ^ 

on the annihilation operators a(0) and b(0). The next 
step then consists in noticing that for <j) = ipt, lu = St, 
and £ = — in*t, the operator 

T 2 (k, 8; t) = e - it [ K * atb M 5 / 2 H ata+btb+1 )+ Kab ] ^4) 

is responsible for generating the terms between braces in 
the solutions a(t) and b(i) - see equations (fTU| and (fiTI) . 
After these considerations, it is easy to show that 

U(t) :=R(w ,w 6 ,<5;t)T 2 M;i) (15) 

provides the solutions obtained above by means of the 
mathematical operation sketched in the identities a(i) = 
Ut(*)a(0)U(t) and b(t) = U t (i)b(0)U(t). Note that the 
particular nonresonant system here studied represents a 
first dynamical application of the results obtained in Sec- 
tion II for the two-mode bosonic representation, where 
the connection between the time-evolution operator (|15p 
and the unitary operator ([9]) is promptly established. It 
is worth emphasizing that different physical systems can 
also be used for explaining the dynamical origin of the 
real parameter oj present inT 2 ((,w) (e.g., see Ref. [13), 
but this fact still deserves be carefully investigated. 

IV. AN INITIAL STUDY ON ENTANGLEMENT 
VIA WIGNER FUNCTION 

Let us initially consider a specific subset of bipartite 
physical systems described by continuous variables such 
that the parametric amplifier model here studied (or the 
simple model for parametric frequency conversion pro- 
posed by Louisell et al. [irj| , and subsequently studied in 
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detail by Tucker and Walls [19|]) constitutes a particular 
element. Furthermore, let us also assume that p(t) = 
U(t)p(0)LF(t) describes the dynamics of any bipartite 
system belonging to this subset whose unitary time- 
evolution operator U(t) is related to the Hamiltonian 
operator H(t); by hypothesis, the initial density operator 
p(0) represents the system prepared at the initial instant 
t — in any disentangled (entangled) state. The symme- 
tric characteristic function for this class of bipartite phy- 
sical systems is given by "^(G; t) := Tr[D(G)p(i)], where 
D(G) := exp(— G^EO) defines a displacement operator 
written in terms of the matrices G* = ( £a £a Q &> ) > 
E = I ® $ with I = diag(l, 1) and S = diag(l, -1) being 
the respective 2x2 identity and symplectic matrices, and 
©t = ( at a b ) . Note that due to the cyclic pro- 
perty of the trace operation, the symmetric characteristic 
function can also be evaluated through the mathematical 
statement Tr[U t (£)D(G)U(i)p(0)], where the time evo- 
lution of the displacement operator plays an important 
role. 

Next, let us suppose that the action of the unitary 
time-evolution operator on the matrix © transforms the 
annihilation (creation) operators a (a.') and b (b^) fol- 
lowing the general rule H (f) = Ut(i)OU(i) = T(i)0, 
where 



T(t) = 



Kit) 



Vb{t) 



Xa(t) 

Xb(t) 

rim 



Va(t) 
X*a(t) 
Tlb(t) 



represents a 4 x 4 matrix whose elements are c-number 
functions determined from the Heisenberg equations for 
© H (i) with specific initial conditions that preserve the 
quantum mechanics (the subscript H indicates that ope- 
rators are in the Heisenberg picture). For instance, when 
t = the c-number functions should necessarily imply in 
the mathematical identity T(0) = diag(l, 1, 1, 1); while 
for t > the commutation relations 

[a H (t),4(*)] = 1= [b H (t),bt(t)] =1, 
[a H (i),b H (t)] = 0, [a H (i),bt(i)] =0, 

provide extra relations for the c-number functions which 
permit us to solve completely the Heisenberg equations, 
namely 

\Va(t)\ 2 -W a (t)\ 2 + \Xa(t)\ 2 -\Va(t)\ 2 = 1, 

\v b (t)\ 2 -\v b (t)\ 2 + \ Xb (t)\ 2 -\v b (t)\ 2 = 1, 

Ha(t)v b (t) - V a (t)Hb{t) + Xa(t)f]b(t) ~ Va(t)Xb(t) = 0, 

M*K(*) - "»(*K(<) + xa(t)xt(t) - Va(tK(t) = o. 

The first immediate consequence of these results is that 
TJt(£)D(G)U(t) will produce a new displacement opera- 
tor D(Y) = exp(-YtEO) with Yt = ( (3* a (3* b b ), 
where the new elements are connected with the old ones 
through the equalities 

Pa = H* a (t)£ a - Va{t)Ca + - V b {t)(,t, 

Pb = X* a (m a -ria(t)C+xt(m-r)b(m, 



plus their respective complex conjugates. For conve- 
nience in our calculations, /3 a t b ) means a short nota- 
tion for /3 a (&) = P a (b) (£a, Cfci t) ■ It is worth noticing that 
tf(G;t) may now be written as ^(Y; 0) = Tr[D(Y)p(0)] 
(i.e., the symmetric characteristic function ^(G; t) is thus 
specified in terms of the form it takes at t — 0, which cor- 
roborates the results obtained by Mollow and Glauber 
[l6| for the parametric amplifier model), and this fact 
will bring some insights into the study of entanglement 
for bipartite systems via Wigner function. 

The Wigner function for this particular subset of bi- 
partite physical systems may then be defined as the four- 
dimensional Fourier transform of the symmetric charac- 
teristic function ^(G;i), that is, 



d 2 £ a d 2 tb 



exp (G t EX)^'(G;i), (16) 



with ID := ( a* a a a I a b ). Note that the variables of 
integration £ a and can be changed to (3 a and (3 b in this 
equation, once the Jacobian matrix of the transformation 
has a determinant whose absolute value is equal to one 
- this fact implies in the identity cP£ d 2 £& = d 2 (3 a d 2 (3 b . 
Hence, the integration can now be conveniently carried 
out through the mathematical relation 



= 3T(Z;0), 



(17) 



where ffl := (7* 7a l b 7b ) defines a new matrix with 
elements given by[30( 

7 a = fi*(t)a a - v a (t)a* + i-L* b {t)a b - v b {t)a* b , 
lb = X* a {i)a a -Va(t)a* a +x b {t)a b -T] b (t)al, 

and their respective complex conjugates. Thus, equa- 
tion (frT|) asserts that W{lL;t) can also be expressed in 
terms of the form it takes at the initial time t = 0, since 
the new variables 7 a (a a , a b ; t) and 7b (a a , a ; t) carry the 
information of the dynamical entanglement between the 
variables a a and ab- Furthermore, this result permits us 
to show that for a given bipartite system initially pre- 
pared in any entangled state, it will remain entangled for 
all t > 0; otherwise, if the density operator is described 
at t = as p(0) = p a (0) <& p b (0) (disentangled state), 
the appearance of entanglement in the Wigner function 
will depend exclusively on the dynamics provided by the 
Hamiltonian operator H(t)[3l[ - here associated with the 
time-dependent matrix Z. Following, let us apply the 
results obtained until now to the parametric amplifier 
model, where different initial states of the electromag- 
netic field will be considered. 

Now we focus our efforts in evaluating #^(Z; 0) for 
a well-known family of two-mode electromagnetic fields 
where, in particular, both modes are initially prepared in 
the coherent states, number states, and thermal states, 
respectively. For instance, the symmetric characteristic 
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functions in these situations are given by 

<^ coh (Y ; o) = e ~^l' 3 »l 2 +l' 3 ''l 2 )+ 2i [ Im (' 3 »C)+im(/3 i) C)] 5 
%(Y;0) = e-Ml^l 2 +IAI 2 )L„ a (|/3 a | 2 )i„ 6 (|/3 fc | 2 ), 
% h (Y;0) = e -i[d+^ a )\0 a \ 2 +(i+2n b m\ 2 ] i 

while their respective Wigner functions can be expressed 
as 

W coh (%;0) = 4e- 2 ( | ^- c " |2+l76 - c '' 12 ), (18) 

W n (Z;0) = 4(-l)"^e- 2 (l^l 2+ l^l 2 )L„ a (4| 7a | 2 ) 

x L nb (4\ lb \ 2 ), (19) 
W th (Z;0) = A[(l + 2n a ){l + 2n b )]- 1 

x e -2[(l+2fi a )- 1 | 7a | 2 +(l+2« il )- 1 |7 6 | 2 ] i / 2Q s 

with L n (z) denoting a Laguerre polynomial and n a ( b ) 
being the mean number of photons for a chaotic light 
field. Note that if we consider the solutions obtained 
from the Heisenberg equations for the parametric am- 
plifier model (see previous section), it is easy to show 
that fJ- a (t), r] a (t), v b {t), and Xb(t) (once the further time- 
dependent coefficients do no exist) determine completely 
the c- numbers /3 a(b) (£ a , £ 6 ; t) and 7 o(b ) (a a , a b ;t). Indeed, 
this last step leads us to characterize precisely "if (Y; 0) 



and W(Zi;0). Moreover, when <5 



0, K G 



and 



(uj a + ui b )t = 3ir/2, Eq. (fTSjl coincides with #ip R (7a, 7&) 
f° r Ca(fc) = (two-mode vacuum state), this function be- 
ing that used by Braunstein and Kimble [12] in the theo- 
retical description of teleportation involving continuous 
quantum variables. [32j| 

Finally, let us say some few words about the results 
obtained in this section. It is important to emphasize 
that equation l|17p permits us to describe qualitatively 
the entanglement of a specific subset of bipartite sys- 
tems, where now we can promptly identify two distinct 
origins of this quantum effect: the first associated with 
the entanglement in the initial conditions (here labelled 
by kinematical entanglement), while the second is res- 
ponsible for the dynamical entanglement via Hamiltonian 
operator. In addition, the normally ordered moments 



a^(t)a'(i)b^(t)b s (i))=r^^e 



(p,q,r,s) i(|5 a | 2 + |? b | 2 ) 



(21) 



with 



•£a,£j,£b,« 



QP+q+r+s 



and {p, q, r, s} € N, can also be used to investigate some 
recent proposals of inseparab ility criteria for continu- 
ous bipartite quantum states [U [23|. In the next sec- 
tion, we will establish a reasonable measure of entangle- 
ment which is based on the results obtained by Piatek 
and Leohski [lil ] for the intermode correlations in phase 
space. 



V. ENTROPY FUNCTION ALS FOR 
CONTINUOUS PHASE-SPACE 
REPRESENTATIONS 

In order to establish a quantitative estimate of entan- 
glement for bipartite systems, let us introduce some basic 
definitions of entropy functionals in continuous phase- 
space representations. The first definition is based on 
the joint entropy fill ] 



E[jT;t] 



CpCtaCpCtb 



Jf{a a , a b ; t) In [Jt?(a a , a b ; t)] 



(22) 

where Jif(a a ,a b ;t) := (a a , a b \p{t)\a a , a b ) denotes the 
Husimi function in the continuous coherent-state repre- 
sentations for a bipartite system described by the density 
operator p(t). Note that Fi[Jff; t] presents certain proper- 
ties inherent to its definition which deserve be mentioned: 
(i) the probability distribution function Jff(a a , a b ; t) is 
strictly positive and limited to the interval [0,1]; con- 
sequently, the joint entropy (|22|) characterizes a well- 
defined function whose behaviour does not present any 
mathematical inconsistencies. Moreover, (ii) this defi- 
nition essentially measures the functional correlation be- 
tween the continuous variables used to describe each part 
of the joint system. Hence, E[Ji?; t] can be considered as 
a natural extension of that definition employed by Wchrl 
[l~2| for information entropy. [33j 

The second definition consists of functionals related to 
the partial entropies 



E[j4? (A) ;t} := - 



d 2 a a 



Jf< A >(a a ;t) In \,^ K \a a -t) 

(23) 



and 



E[^f (B) ;<] := - 



d 2 a b 



^ (B) K;t)ln ,yf {B} (a b ;t) 



which depend basically on the marginal Husimi functions 



^ {B) (a b ;t) = 



d 2 a b 



d 2 a a 



J4?(a a ,a b ;t), 



Jf(a a ,a b ;t). 



Such marginal Husimi functions carry information of the 
entanglement between the subsystems 'A' and 'B', since 
the partial trace over the continuous variables a a ^ of 
a particular subsystem 'A' ('B') allows the introduction, 
via time-evolution operator U(i) and/or initial density 
operator p(0), of important correlations in the bipartite 
states. 

Let us derive now some mathematical relations among 
these entropy functionals from the Araki-Lieb inequality 
El, i.e., 



E[J>f (A) ;£] -E[Jf (B) ;t] 



< E[JT;i] < Epf (A) ;i] 
+ E[JfT^;t]. (25) 
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For instance, the rhs of this inequality corresponds to the 
subadditivity property for the Wehrl's entropy functio- 
nals, while the equal sign reflects the complete disentan- 
glement between the subsystems 'A' and 'B' of the joint 
system. In addition, the conditional entropies (Til . Il2| 

E[JT / Jf (A) ; t] = E[Jt; t] - E[JT (a) ; t] (26) 

and 

E[Jf / Jf < B > ; t] = E[JT; t] - E[Jf ^ ; t] (27) 
lead us not only to establish the balance equation 

E[Jf/J^ (A) ;t}+E[J^ {A) ;t} = E[J^/jT (b) ; t]+E[Jf (B) ; t], 

but also to determine the further inequalities 
E[Jf/Jif M ;t} < E[Jf (B) ;i] 

and 

E[Jf/jT (B) ;i] < E[Jf (A) ;t] 

from the subadditivity property. It is worth noticing 
that the equal signs hold in both situations only when 
the c-numbers a a and a b are functionally uncorrelated, 
namely, the bipartite Husimi function J^(a a , a b ;t) fac- 
torizes in the product of the marginal Husimi functions 
,^ A \a a -t),^\a b -t). 

In what concerns the class of entropic functionals listed 
until now, it is convenient to define a new functional li- 
mited to the closed interval [0, 1] for any t > 0, which 
permits us to avoid any ambiguity in the significance of 
the subadditivity property. Thus, let us introduce the 
correlation functional 

(29 ' 

which can be used to measure the functional correlation 
(that is, the 'degree of entanglement') between the parts 
'A' and 'B' of the joint system, having as reference a fac- 
torizable bipartite Husimi function J^(a a , a b ) t). In this 
definition, the inferior limit of the closed interval [0, 1] 
corresponds to the situation where the bipartite Husimi 
function factorizes, while its superior limit reflects the op- 
posite situation (i.e., when the bipartite states are maxi- 
mally correlated); besides, the global factor '2' is associ- 
ated with the number of possible combinations between 
the constituent parts of the system under investigation. 
As a first practical application of these entropy functio- 
nals we will consider two different examples of initially 
uncoupled bipartite states for the parametric amplifier 
model. 



A. Coherent states 



t = of the model governed by the Hamiltonian (|12p . 
Here, |Ca,C&) = ICo) ® ICfo) where |C a (6)) characterizes the 
coherent states related to the signal (idler) mode of the 
electromagnetic field. Thus, after some nontrivial alge- 
bra, the bipartite Husimi function assumes the closed- 
form 

x e 2H4A;(7.-C.)(T*-6)] (30) 

such that 

7«(b) = e l{ ^ +s/2)t [cosh(pi) - i{6/2tp) sinh(^)] a a(b) 
+ i(«7 V )e- i ^w+«/a)t Binh( ^ )a j (a)j 

^ i(ft*/V) sinh(yrt) 

cosh(^) + i(8/2ip) sinh(^t) ' 

A a = [cosh(v3i) +i(S/2(p)smh((pt)]~ 2 . 
Furthermore, the marginal Husimi functions 

■K^(a a ;t) = \A \e-\ A ^-^ 2 (31) 

and 

^(a b ;t) = \A \e-^ a ^ 2 , (32) 

with 

e ( 6 ) = | [cosh(^t) + i{5/2ip) sinh(^)] Ca(b) 

-i( K 7^)sinh(^)C (Q) } e - i ^< & » +5/2) *, 

permit us to verify how the dynamical entanglement in- 
troduces correlations between the continuous variables 
Ca(b) an d Cb(a) through the label c a (py Indeed, for t = it 
is possible to demonstrate that J4? co h{a a , a b ; 0) factorizes 

in the product J l fj^(a a ;0)^^(ab;0)- However, if one 
considers t > 0, the dynamics governed by the process 
of parametric amplification introduces significant corre- 
lations between the signal and idler modes. Following, let 
us quantify these correlations with the help of Eq. ([29)) . 

The Gaussian structures present in the bipartite and 
marginal Husimi functions allow us to obtain closed form 
expressions for the respective joint and partial entropies, 
namely, 

E[Jf coh ;t] = 2 + In [1 + (|k|/</?) 2 sinh 2 (^)] (33) 

and 

E[^ c ( h B) ; t] = 1 + In [1 + (M/v?) 2 sinhV)] . (34) 

Consequently, the correlation functional can be promptly 
evaluated as follows: 



In this first example, we consider the bipartite pure 
states p(0) = |Ca, (b){(a, £b\ as being the initial state at 
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Besides, if one considers the limit ipt 3> 1 for sufficiently 
long time, the correlation functional goes asymptotically 
to 1, this value being associated with the specific dyna- 
mics given by the Hamiltonian fT2")) . On the other hand, 
for (pt = we obtain C[^ h; 0] = 0, which corroborates 
the factorization process of the bipartite Husimi function 
J^oh(a a , a b ;0). 



B. Thermal states 

Now, let us suppose that both signal and idler modes 
were initially prepared in the thermal states. In this case, 
the initial density operator admits the expansion 

n a n b J 

where n a i b ) represents the average occupation number of 
each mode, and P(£ , ( b ) = \( a , C&)(Ca, Cb| is the projector 
of the coherent states. This typical example of bipartite 
mixed states constitutes an important starting point in 
the investigation process of the dynamical entanglement 
due to the parametric amplifier model. So, after some 
calculations, the bipartite Husimi function can be written 
as 

J% h (a a ,a b ;t) = [0^(1, l;*)]" 1 e -(^l 2 +^l 2 ) 

a 2Ke(£ ai ot a a b ) /gg\ 



x e 



with 



.,(0,1;*) 



^n a ,n b (1? 1; *) ^n a ,n b ( 1 ; 1 5 ^) 

l ab = i(«;/V) smh(<pt) [cosh(ipt) — i(6/2(f) smh(ipt)] 
x (n„ + n b + 1) [5^(1, 1; t)}' 1 j^+^+W, 

and 

^n a ,n b (x,y;t) = (n a x + l)(n b y + 1) 

+ (n a + n h + 1)(|k|/</?) 2 sinh 2 (^); 

while its joint entropy assumes the closed-form 

E[J«| h ;t]=2 + ln[Sfn. lfii (l ) l;t)]. (37) 

In addition, the marginal Husimi functions 

and 

^ B) (a h ;t) = [^ 



l l~ |2 

(38) 



permit also to determine explicitly the partial entropies 



E[^ A) ;i] = 1 + ^^,^(1,0;*)], (40) 
E[^ B) ;t] = l+ln[Sf« a , fl6 (0,l;t)]. (41) 



In this moment becomes important noticing that the sub- 
additivity property is not violated for any t > 0, the equa- 
lity E[^ h ; 0] = E[J^ A) ;0] + E[J^ b) ;0] being consistent 
with the factorization of J^th (a , a b ; 0) in the product 



^ A) (a a ;0)^(a 6 ;0). 

Following, let us investigate the correlation functional 
C[J£th; t] for * > 0, since C{J#th', 0] = reflects the disen- 
tanglement between the bipartite mixed states. In fact, 
for t =/= this functional assumes any values into the open 
interval (0, 1), the maximum value 1 being reached only 
for sufficiently long time. Hence, we can conclude that: 

(i) any initially uncoupled states become entangled in the 
process described by the Hamiltonian (jT2j) , which con- 
firms the results obtained by Dodonov et al. [2l| through 
the evaluation of the inverse negativity coefficient (this 
measure also leads us to estimate quantitatively the 'de- 
gree of entanglement' of a particular class of bipartite 
states such as those studied in this section); furthermore, 

(ii) the dynamics here characterized by the physical sys- 
tem under investigation produces a maximum estimate 
of entanglement for any initially disentangled states, this 
value (namely, 1) being considered as a specific quantum 
signature of the parametric amplifier model. 



VI. CONCLUSIONS 

We have established a set of interesting formal results 
within the scope of quantum optics and quantum infor- 
mation theory that allows us, among other things, 

• to define a new class of unitary squeezing transfor- 
mations related to su(l, 1) Lie algebra which gene- 
ralizes, for certain particular representations of its 
generators, the one- and two-mode squeezing ope- 
rators 3; 

• to discuss, from the physical point of view, how 
the generalized two-mode squeezing operator can 
be generated through a slightly modified version of 
the parametric amplifier model [Tol ]; 

• to obtain a general integral representation for the 
bipartite Wigner function whose integrand is ex- 
pressed as a product of two terms which are res- 
ponsible for the dynamical and kinematical entan- 
glement; and finally, 

• to estimate quantitatively the 'degree of entangle- 
ment' related to an ideal bipartite system (we are 
discarding the unvoidable coupling with the envi- 
ronment in this context) by means of a theoretical 
framework [ll[ which is based on the Wehrl's ap- 
proach fl2j for the entropy functionals. 



In fact, equation (|29|) and its respective properties men- 
tioned properly in the body of the text permit us to es- 
timate, through an entropic approach, the entanglement 
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effects for a wide class of electromagnetic field states, in- 
cluding Gaussian and non- Gaussian states. As a con- 
cluding remark, it is worth mentioning that these re- 
sults have also potential applications in modern research 
on quantum teleportation [22|, quantum tomography 
0, 0] , and quantum computation [27| (once continuous- 
variable entanglement can also be efficiently produced 
using squeezed light and linear optics 28]), as well as on 
the foundations of quantum mechanics through its exten- 



sions to the su(2) Lie algebra [2S 
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